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Abstract 

It has been argued that the emergence of roughly periodic orientation preference maps (OPMs) in the primary visual cortex 
(VI) of carnivores and primates can be explained by a so-called statistical connectivity model. This model assumes that input 
to VI neurons is dominated by feed-forward projections originating from a small set of retinal ganglion cells (RGCs). The 
typical spacing between adjacent cortical orientation columns preferring the same orientation then arises via Moire- 
Interference between hexagonal ON/OFF RGC mosaics. While this Moire-Interference critically depends on long-range 
hexagonal order within the RGC mosaics, a recent statistical analysis of RGC receptive field positions found no evidence for 
such long-range positional order. Hexagonal order may be only one of several ways to obtain spatially repetitive OPMs in 
the statistical connectivity model. Here, we investigate a more general requirement on the spatial structure of RGC mosaics 
that can seed the emergence of spatially repetitive cortical OPMs, namely that angular correlations between so-called RGC 
dipoles exhibit a spatial structure similar to that of OPM autocorrelation functions. Both in cat beta cell mosaics as well as 
primate parasol receptive field mosaics we find that RGC dipole angles are spatially uncorrelated. To help assess the level of 
these correlations, we introduce a novel point process that generates mosaics with realistic nearest neighbor statistics and a 
tunable degree of spatial correlations of dipole angles. Using this process, we show that given the size of available data sets, 
the presence of even weak angular correlations in the data is very unlikely. We conclude that the layout of ON/OFF ganglion 
cell mosaics lacks the spatial structure necessary to seed iso-orientation domains in the primary visual cortex. 
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Introduction 

Many neurons in the primary visual cortex (VI) respond 
preferentially to edge-like stimuli of a particular orientation [1]. In 
carnivores and primates, orientation preference exhibits a 
columnar arrangement such that neurons positioned on top of 
each other from the white matter to the pia typically prefer similar 
orientations. Tangential to the visual cortical layers, orientation 
preference changes smoothly and progressively [1] except at the 
centers of so-called prnwheels where neurons exhibiting the whole 
range of orientation preferences are located in close vicinity [1,2]. 
The progression of orientation preferences across the visual 
cortical surface (Orientation preference map, OPM) appears as 
organized by a semiregularly spaced system of pinwheels and 
adjacent columns preferring the same orientation over roughly a 
millimeter distance [3-11]. 

Most models for the emergence of OPMs during postnatal 
development assume that their layout is determined by intracor- 
tical mechanisms (e.g. [12-17]). However, several recent studies 
advance the notion that the structure of OPMs may result from a 
statistical wiring of feed-forward inputs from the mosaic of ON/ 



OFF retinal ganglion cells (RGCs) to VI [18-21] (Fig. 1A), an idea 
pioneered by Soodak [22,23]. ON/OFF ganglion cells are 
arranged in semiregular mosaics across the retina and project to 
the lateral geniculate nucleus (LGN) of the thalamus. Thalamic 
receptive fields resemble RGC receptive fields in shape, size, and 
spatial distribution [24,25]. The retinotopic map [26-29] allows 
neighboring retinal/ thalamic ON and OFF center cells to project 
to neighboring neurons in the primary visual cortex. Most nearest 
neighbor RGCs are ON/OFF pairs [30]. According to the 
statistical connectivity model, a VI neuron predominantly samples 
feed-forward inputs from geniculate projections in its immediate 
vicinity [31]. If so, it is likely to receive input from a single pair of 
ON/OFF RGCs, a so-called dipole [20] (Fig. 1A, left). The 
neuron's receptive field would then be dominated by one ON and 
one OFF subregion (Fig. 1A, middle left) and its response 
orientation-tuned [18-20]. In this picture, the preferred orienta- 
tion is determined by the orientation of the RGC dipole (Fig. 1 A, 
middle right, right). Consequently, a key prediction of the 
statistical connectivity model is that orientation preference across 
the surface of the primary visual cortex should mirror the spatial 
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Figure 1 . RGC dipoles and the statistical wiring model according to [20]. A Most left: A dipole of an ON center (empty circle) and OFF center 
(filled circle) retinal ganglion cell (RGC). The black line connecting the two cells indicates that the two cells form a dipole. A VI cell with input 
dominated by this dipole has a receptive field with side-by-side subregions of opposite sign (middle left) and is consequently tuned to a specific 
orientation (j> (middle right). We represent the preferred orientation of the V1 cell by the color of the bar connecting the two RGCs (most right). Note 
that the preferred orientation of the V1 is orthogonal to the bar connecting the two RGCs. B The statistical connectivity model for orientation 
preference maps. The receptive field midpoints of ON/OFF center RGCs are arranged in semiregular mosaics. The input to a cortical cell is dominated 
by a single pair of ON/OFF dipole and the cortical units have oriented receptive fields. If RGC dipole orientations are locally correlated, orientation 
preference 0 within layer 4 of V1 is predicted to vary smoothly resulting in a smooth and continuous map of orientation preferences. C Parametrized 
definition of RGC dipoles. ON/OFF pairs with distance smaller than a parameter d are considered dipoles (black lines). For the centered OFF cell, 
preferred orientations of dipoles are indicated. With increasing d, the number of dipoles increases and one RGC can form multiple dipoles. 
doi:10.1371/journal.pone.0086139.g001 



distribution of the ON/OFF dipole angles in the RGC mosaics 
[18,20,22,23] (Fig. IB). 

Paik and Ringach [20,21] showed how this model can explain 
the experimentally observed roughly periodic structure of visual 
cortical OPMs. In their model, ON and OFF cell mosaics are 
assumed to form two independent noisy hexagonal lattices. 
Superposing these two lattices leads to a hexagonal pattern of 
dipole orientations via Moire interference [32]. The statistical 
connectivity hypothesis then implies that this hexagonal pattern is 
mapped onto the cortex creating a roughly hexagonal OPM [20]. 

Hore et al. recently cast substantial doubts on the presence of 
hexagonal order in experimentally measured RGC receptive field 



mosaics [33]. They found that the positions of receptive fields of 
ON/OFF RGCs in monkey retina are inconsistent with long- 
range hexagonal order and are much better described by a so- 
called pairwise interacting point processes (PIPP) in a parameter 
regime where long-range positional correlations are absent [33— 
35]. Such PIPP mosaics lack the long-range order necessary to 
create a Moire-Interference pattern and hence OPMs predicted by 
the statistical connectivity model with PIPP mosaics do not exhibit 
the experimentally observed spatially repetitive arrangement of 
orientation columns [18,19,33]. 

Moire-Interference of hexagonal RGC lattices constitutes one 
particular way of creating an ordered arrangement of regularly 
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spaced orientation columns in the statistical connectivity model. 
Other spatial arrangements of RGCs are conceivable that might 
lead to spatially repetitive cortical orientation maps resembling the 
ones experimentally measured. Therefore, the lack of hexagonal 
structure in RGC mosaics found by Hore et al. does not per se 
dismiss the statistical connectivity hypothesis. 

Here, we investigate a fundamental requirement on the spatial 
structure of RGC mosaics to seed the emergence of spatially 
repetitive cortical OPMs: a spatial correlation of RGC dipole 
angles across the retina. RGC dipole angle correlations are 
predicted to exhibit a spatial structure similar to that of OPM 
autocorrelation functions. This means that RGC dipole angles 
have to be locally positively correlated and anti-correlated on 
intermediate scales. The precise values for both of these scales 
depend on the column spacing of the OPM as well as the cortical 
magnification factor. We first systematically analyze two previ- 
ously published mosaics of cat beta cell somata positions [30,35] as 
well as one primate parasol mosaic of RGC receptive field centers 
[36] with respect to their dipole angle correlation functions. In 
both species, we are unable to detect any statistically significant 
positive or negative correlation. Since all three mosaics analyzed 
each contain only around 100 cell positions (or RFs center 
positions), the absence of detectable correlations might be a 
consequence of the small size of the data sets. To address this 
question, we introduce a novel point process that generates 
mosaics with a tunable degree of spatial correlations of dipole 
angles. The spatial structure of these angular correlations is 
designed such as to match the autocorrelation functions of 
experimentally measured OPMs. At the same time, the local 
spatial statistical properties of the resulting model RGC mosaics 
agree well with the experiment. On the one hand, the mosaics 
generated by this process demonstrate that hexagonally organized 
RGCs are indeed not necessary to obtain realistic OPMs within 
the statistical connectivity model. On the other hand, the mosaics 
generate by the point process can be used as reference mosaics to 
access the amount of data needed to detect the presence or 
absence of correlations. Finally, by statistical comparison of 
mPIPP model mosaics and data, we show that, given the size of 
our data sets, the presence of even weak angular correlations in the 
data can be ruled out. We conclude, that the layout of ON/ OFF 
ganglion cell mosaics apparently lack a fundamental feature 
necessary to explain the emergence of spatially repetitive 
orientation preference maps in VI within the dipole approxima- 
tion of the statistical connectivity model. Our results suggest that 
the ordered arrangement of orientation columns is unlikely to 
originate from the spatial layout of RGCs and, hence, likely set by 
intracortical mechanisms. 

Results 

Dipole Orientation Correlation Function in Cat and 
Primate RGC Mosaics 

We first determined the correlation function of dipole orienta- 
tions in the two published cat beta cell mosaics. These mosaic 
fields will be referred to by their keys: w81sl and m623. Field 
w81sl was created by digitizing the map shown in Figure 6 of [30]. 
Field m623 was taken from [35]. To limit the restriction of 
considering only nearest neighbor ON/OFF cells pairs, we 
followed the flexible dipole definition introduced in [37] 
(Fig. 1C). In this definition, a parameter d is introduced, 
describing a distance below which neighboring pairs of ON/ 
OFF cells are considered as dipoles (see Eq. (4)). The larger d, the 
more dipoles are formed by each RGC cell (cf. Fig. 1C). The 
nearest neighbor distance distribution of the mosaic defines a 



range of sensible rf-values for each mosaic. For instance, the 
nearest neighbor distribution of cat beta cell mosaics peaks around 
60fim [30] and, therefore, rf-values smaller than 60fim lead to the 
extraction of only very few dipoles. On the other hand, values 
larger than 100 fim lead to many dipoles between non-nearest 
neighbors. 

Figures 2A-C show the m623 mosaic along with all the dipole 
pairs, color-coded according to their orientation extracted for 
d = 60 pirn (A), d = 80fim (B), and d = lOOfim (C) (see Materials and 
Methods). Figures 2D-F display the dipoles found for the mosaic 
w81sl. For d = 80fim and 100/ra, dipoles in both mosaics appear 
as organized into linear chains which, at first sight, one might take 
as an indication of spatial correlation in the dipole orientations 
[37]. However, a quantitative analysis reveals that such correla- 
tions are absent for all d-values (Figs. 2G-J). The only statistically 
significant correlation present in both mosaics is a weak anti- 
correlation on very short scales (<100/ra) for d=W0[im. This 
finding will be explained in more detail below. 

The two mosaics analyzed so far are based on the position of 
cell bodies and not those of the actual RGC receptive fields. Since 
these two do not match necessarily, it is important to test whether 
RGC receptive field mosaics substantially differ from cell position 
mosaics with respect to their dipole correlation structure. 
Therefore, we next repeated the above analysis for a previously 
published primate parasol cell receptive field mosaic (referred to 
by its key G09 [36] in the following). Figures 3A-C show the G09 
mosaic along with all the dipole pairs, color-coded according to 
their orientation extracted for d = 60 fim (A), d = %QpLm (B), and 
d= 100 fim (C). Again, a quantitative analysis reveals that dipole 
angular correlations are absent for all d-values (Figs. 3D-F). 

To analyze local correlations more systematically, we varied the 
parameter d for all three mosaics in the data set over the entire 
range of sensible values and determined angular correlations in 
each of the mosaics in the first distance bin (see Materials & 
Methods). Figures 4A and B depict the results of this analysis. 
Interestingly, for small rf-values, all mosaics exhibit a weakly 
positive local correlation, whereas for larger d-values dipoles 
appear anti-correlated in all three mosaics (see Fig. 2J for an 
example). Intuitively, this dependency of correlation values can be 
understood is a consequence of our flexible dipole definition 
together with a typical distance between nearest neighbor RGC 
(Fig. 4A). For small and intermediate <i-values, an ON-cell 
positioned between several OFF-cells forms dipoles with mosdy 
one or two of them. If dipoles are formed with two OFF cells, due 
to the regular spacing of OFF cells, the angles formed by these two 
dipoles are likely positively correlated (Fig. 4A, left inset). For 
larger af-values, more than two dipoles are typically formed. In this 
case, the regular spacing leads to an effective anti-correlation 
between their angles (Fig. 4A, right inset). To investigate whether 
such correlations suffice to seed the development of spatially 
repetitive and smooth OPMs in V 1 , we compared the correlation 
traces found in the experimentally measured mosaics to model 
mosaics obtained by a pairwise interacting point process (PIPP) 
[33,34]. The PIPP is a method for the generation of a spatial 
distribution of points specifying only pairwise interaction between 
individual points. It has been previously shown to accurately 
reproduce the spatial statistics of experimentally measured RGC 
mosaics [33,34]. With parameters fitted to experimental data, the 
PIPP generates regularly spaced RGC mosaics with radially 
isotropic autocorrelograms and lack of long-range positional order 
[33]. For such mosaics, the statistical connectivity framework 
predicts OPMs that lack a typical column spacing and are 
qualitatively different from experiment [18,33]. Figures 4D,E 
show that the local correlation of dipole angles in the model 
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Figure 2. Spatial correlations of dipole orientations are absent in cat beta cell mosaics. A ON/OFF cells (empty/filled circles) for the cat 
beta cell mosaic m623 [35]. Preferred orientation of dipoles extracted for d = 60/Pn are shown as colored bars. Colorcode as in Fig. 1 C. B as A but for 
rf = 80/»j). C as A but for d = YOOfUn. D-F as A-C but for cat beta cell mosaic w81s1 [35], G Correlation of dipole orientations for m623 (red) and w81s1 
(blue), calculated from dipoles extracted for rf = 60/»;i. Error bars indicate 95% confidence intervals of bootstrap distributions. H as G but for 
d=i0/im. J as G but for d=\00^im. 
doi:1 0.1 371 /journal.pone.00861 39.g002 



mosaics exhibit the same d-dependence as the data. For small d- 
values mosaics exhibit a weakly positive local correlation, whereas 
for larger rf-values dipoles are weakly anti-correlated. We conclude 
that the weak local correlations found in experimental mosaics are 
a consequence of our flexible definition of RGC dipoles together 
with a typical spacing between neighboring cells. The quantitative 
match between PIPP mosaics and experimental data together with 
the fact that PIPP mosaics seed unrealistic cortical OPMs in the 
statistical connectivity model indicate that these correlations are 
not sufficient to explain the layout of realistic cortical OPMs. 



An OPM-modulated Pairwise Interacting Point Process 

One might wonder whether the absence of detectable positive 
or negative correlations in dipole angles is merely a consequence of 
the small sizes of each of these data sets. In fact, all three mosaics 
analyzed each contain only about 100 cell positions (or RF center 
positions) and a similar number of dipoles. To clarify whether such 
small data sets are sufficient to detect both relevant dipole 
correlations, it is necessary to design model mosaics, e.g. denned 
by a some point process, with realistic spatial statistics and with a 
known degree of dipole angle correlations. Using such a point 
process, one can then ask whether the size of the available data sets 
in principle permits detection of such correlations. However, up to 
this point, a suitable point process with a known degree of dipole 
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Figure 3. Spatial correlations of dipole orientations are absent in a primate parasol receptive field mosaic. A ON/OFF cells (empty/ 
filled circles) for primate parasol cell receptive field mosaic G09 [36]. Preferred orientation of dipoles extracted for rf = 60/im are shown as colored 
bars. Colorcode as in Fig. 1C. B as A but for <:/ = 80/im. C as A but for d = lOO^im. D Correlation of dipole orientations for mosaic G09, calculated from 
dipoles extracted for <f = 60/(m. Error bars indicate 95% confidence intervals of bootstrap distributions. E as D but for d =80/im. F as D but for 
d = lOOjum. 

doi:1 0.1 371 /journal.pone.00861 39.g003 



angle correlations has not been proposed. We now introduce a 
pairwise interacting point process (PIPP) with such characteristics, 
and start by briefly outlining the definition of the conventional 
PIPP for RGC mosaics as introduced in [34] . 

The PIPP is a method for the generation of a spatial distribution of 
points specifying only pairwise interaction between individual points. 
Interactions between points are usually specified in pairwise 
interaction functions. The product of these pairwise interaction 
functions for a specific location for all possible pairs of points gives 
the probability of finding a point at a particular position. For 
bivariate data such as the positions Xq N and Xq FF of a mosaic of 
ON-cells and a mosaic of OFF-cells, the PIPP is characterized by two 
intra-mosaic interaction functions ^on,on(I x on — x onD> 
/'off.off(I x (3ff ~ x offD an d one inter-mosaic interaction 
' ! on,off(I x on — x offD- To describe the positioning of beta cells in 
the cat retina, Eglen et al. [34] used a parametric form of repulsion 



f 0 if r<8 

h{r)= { 1 ( , r ~ S A T « 
1^ 1 - expl — I 1 J if r>5 



for all three interaction functions, with r = |Xq N — Xq N |, 

>"=I x OFF -X OFfI 0r r= l X ON _X OFFl for A ON,ON, &OFF.OFF Or 

/'on.off, respectively. By fitting the parameters a and q> to 
experimental data, they showed that inter-mosaic interactions are 
sufficiently described by solely ensuring that two cells are not less 
than the soma distance apart, i.e. 

f0 ifr<<5 
hoK,OFF{r) = < . . 1) 

t 1 if r>6, 

with S being the soma diameter. Intra-mosaic interactions were best 
fit by values for a and (p which ensure a semiregular placement of 
RGC cells without long-range positional order [34]. Such mosaics 
lack non-local spatial order in their dipole angles and, hence, the 
statistical connectivity framework predicts orientation maps that lack 
a typical column spacing [18,33]. 

To introduce correlations of dipole angles into the PIPP, we 
start by formalizing the illustration in Fig. 1A. The dipole vector 
Xq N — Xq FF between an ON/OFF pair of RGCs points in the 
direction arg (xJ-, N — Xq FF ) in the interval [— The preferred 
orientation of the dipole (Fig. 1A, most right) can then be 
mathematically expressed as 
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Figure 4. Measured local correlation values depend on choice of dipole distance parameter d. A Local correlations in cat beta cell 
mosaics (red: m623, blue: w81s1) as a function of the dipole extraction parameter d (see Fig. 1C). B As A but for primate parasol cell receptive field 
mosaic G09. C Weak local positive or negative correlations emerge in the same RGC configuration, depending on dipole extraction parameter d. D As 
A but for simulated PIPP mosaics with parameters fitted to cat beta cell mosaics m623 (red) and w81s1 (blue). E As B but for simulated PIPP mosaics 
with parameters fitted to primate parasol receptive field mosaics G09. All error bars indicate 95% confidence intervals of bootstrap distributions. 
doi:10.1371/journal.pone.0086139.g004 



<Mxo N ,x' OFF ) = mod (arg(x^ N - x' OFF ) + n/2,n) , (2) 

varying in the interval [0,7t). The main idea is now to modify the 
inter-mosaic interaction function (Eq. (1)) to add a dipole 
correlational structure that matches the spatial correlations of 
orientation preferences in visual cortical OPMs. A visual cortical 
OPM can be represented as 0(x) with 0e[O,7t). Note that x here 
describes positions on the retina. An OPM measured with optical 
imaging of intrinsic signals [3,38] is naturally given as 0'(K), where 
X specifies cortical position. The retinotopic map x = R(X) 
associates a cortical position X with a position x on the retina. 
Via the inverse transformation, the representation of the OPM in 
retinal coordinates is obtained, i.e. 9(x) = 0'(R ~'(x)). In the 
following, only small subregions of OPMs were considered and the 
retinotopic map was assumed to be linear. 

Using the OPM representation 0(x), we modify the inter-mosaic 
interaction function Eq. (1) by multiplying it with a function 
/! } ,(Xq N ,Xq FF ) that depends on the difference between the 
preferred dipole angle <Kxc, n ,Xq FF ) and the preferred orientation 
specified in the OPM at position (Xq N + Xq FF )/2 (half way 
between the two RGCs): 



-ffoN.OFF (Xqn > X OFF ) = ^ON.OFF (?) 'fly (Xq N ,Xq FF ) . (3) 

We choose 

' I 7^ X ON' X OFf) 

= f exp{ 7 [cos(«xi )N ,x OFF )-«((xi 3N +x{ 3FF )/2)) -l] } if r<d (4) 
\\ if r>d. 

In addition to ensuring that two cells are not less than the soma 
distance apart, this new inter-mosaic interaction function enforces 
that the positioning of an ON cell at position Xq N and an OFF cell 
at position Xq FF is such that the preferred orientation 
</>(Xq N ,x' off ) of the ON/OFF pair is similar to preferred 
orientation specified in the OPM at the corresponding point. In 
this way, the OPM 6(x) is expected to modulate the positioning of 
ON/ OFF cells when numerically simulating the positioning of the 
RGC mosaics with a Monte-Carlo procedure [34] (Fig. ID) such 
that preferred orientations of dipoles align with preferred 
orientations given in the OPM. Note that the preferred orientation 
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of a dipole is orthogonal to the dipole vector Xq N — Xq FF (see 
Fig. 1A). With respect to the PIPP, this model has two additional 
parameters. The strength of modulation through the OPM is 
specified by a modulation parameter y. If y is zero, the positioning 
of ON/OFF pairs is not influenced by the OPM region and the 
model is equivalent to the PIPP model. The larger the value of y, 
the stronger the penalty for OPM and dipole angle differing. We 
refer to the process specified by this inter-mosaic interaction 
function as modulated PIPP (mPIPP). The second parameter is the 
distance d below which neighboring pairs of ON/ OFF cells are 
considered as dipoles (see Eq. (4)). This parameter is taken to be 
the same as the parameter for defining dipoles (cf. Fig. 1C). Again, 
the larger d, the more dipoles each RGC cell is assumed to form 
with surrounding cells (cf. Fig. 1C). The nearest neighbor distance 
distributions of the different mosaics imply sensible values for the 
choice of d. 

We would like to emphasize that by defining the mPIPP as 
above, we by no means imply any influence of cortical orientation 
preference upon the positions of ON/ OFF RGCs during postnatal 
development. The mPIPP merely is a phenomenological algorithm 
to attempt to "reverse engineer" one plausible realization of an 
RGC mosaic with the necessary spatial structure to yield an OPM 
with realistic spatial properties within the statistical connectivity 
model as considered in [20]. 

Statistical Characterization of mPIPP Mosaics 

To characterize the mPIPP defined above, we first investigated 
how the positioning of ON/ OFF RGC dipoles is altered by the 
modulation through realistic OPMs. To this end, we extracted 
small regions of previously published OPMs from cat area 1 7 (data 
courtesy of Z. Kisvarday, see Materials and Methods) and used 
them to modulate the PIPP according to the above description. 
Rectangular OPM regions must be chosen such that the size of the 
retinal region corresponding to the OPM fits the two cat beta cell 
somata position mosaics w81sl and m623 [30,35]. Thus, the size 
of the extracted cortical region depends on the cortical magnifi- 
cation factor mm c /rrirri r at the position the mosaic has been 
recorded from. The larger the magnification factor, the larger the 
extracted cortical region has to be. 

The center of mosaic w8 1 s 1 was located 1 9 below the mid line 
of the visual streak, 4mm from the area centralis [35]. The mosaic 
m623 has been obtained from 5 mm eccentricity and 5.5 below 
the mid line of the visual streak. Both fields are situated at positions 
in the visual field with similar cortical magnification factor (see 
Materials and Methods), estimated to be £= 1.7mm c /mm,., where 
mm c is mm on the cortex and mm, is millimeter on the retina. 
Both RGC fields are of similar linear extent («1 mm). Hence, 
extracted visual cortical regions were of x 1 . 7 mm linear extent 
corresponding to an area of approximately 3 hypercolumns [9,39]. 

Figure 5A show mosaics obtained by simulating an mPIPP with 
a Monte-Carlo procedure for a cortical magnification of 
£= 1.7mm c /mm r , modulated by a cat OPM region for y = 0 
(left), y = 2.5 (middle) and y = 30 (right). ON/OFF pairs less than 
d = 80/jot apart were considered to form dipoles. All other 
parameters were chosen as in [34] for m623 (see also Table SI). 
In the case y = 0 (original PIPP), the dipole orientations are not 
correlated with the preferred orientations specified in the OPM. 
For y = 2.5, some dipoles tend to locally align their orientation to 
match the orientation preference given in the OPM but most 
dipoles have random orientations. For y = 30, the dipole orienta- 
tions are strongly correlated with preferred angles specified in the 
OPM. mPIPP realizations generally display fewer dipoles than the 
experimental mosaic m623 for the same rf-values (cf. Fig. 2B). This 
is true even for the conventional PIPP process (m623: 109 dipoles; 



mPIPP y = 0: 98.6 dipoles on average; y = 2.5: 77.8 dipoles on 
average; y = 30: 55.4 dipoles on average). This can be attributed to 
a slighdy increased inhomogeneity of both the PIPP and mPIPP 
compared to the experimental data (see also Fig. 2 in [34]). 
Figure 5B shows the spatial correlation of dipole angles for all 
three y-values (see Materials and Methods). With increasing y, 
dipole orientations become locally correlated for distance smaller 
than 200/ra. Correlation drops to zero around 200/.i»j distance 
and negative values are obtained between 200 and 500/.ito. 
Moreover, the correlation function of dipole orientations ap- 
proaches the correlation function of the OPM for increasing y 
(solid dark red lines in Fig. 5B). Similarly, the cross-correlation 
between the preferred orientation given by the OPM and the 
orientation of the dipole rapidly increases with y (Fig. 5C). Finally, 
we considered mPIPP mosaics modulated by larger OPM regions 
(size 3x3 column spacings A, Fig. 5D). Again, spatial correlations 
of preferred dipole orientations approach the correlation of 
preferred orientations in the OPM for increasing y (Fig. 5E). 

The dipole angle correlation function periodically modulates 
around zero because in mPIPP mosaics dipole angles roughly 
repeat within a typical distance. This shows that the mPIPP can 
generate RGC mosaics that within the statistical connectivity 
framework predict orientation maps with a typical distance 
between adjacent column preferring the same orientation. The 
typical distance upon which dipole angles roughly repeat depends 
on the column spacing of the OPM used to modulate the PIPP as 
well as the assumed cortical magnification factor. Smaller column 
spacing or larger magnification factors will lead to a smaller scale 
of periodicity in the dipoles and vice versa. Importandy, however, 
dipole angle correlations decay to zero for distances larger than 2— 
3 A, even for large y (Fig. 5E). Hence, long-range spatial order in 
the dipole angles are absent in the mPIPP mosaics. This clearly 
distinguishes the mPIPP dipoles from a pattern of dipoles obtained 
by Moire interference of two noisy hexagonal lattices where long- 
range order is expected [33]. 

We next wondered whether the short-range spatial statistics of 
ON/ OFF cell positions was affected by placing dipoles such that 
their preferred orientation tends to align with realistic OPMs. To 
answer this question, we compared the nearest neighbor distance 
distributions, autocorrelograms of cell positions as well as spatial 
regularity of mPIPP mosaics with the data using statistical 
measures defined in [33,34] (see Materials and Methods). 
Figure 6 depicts the results of this analysis for y = 0;2.5;30 for 
the m623 OFF cell mosaic. Voronoi polygons and autocorrelo- 
grams of all simulated mPIPP mosaics bear close resemblance to 
the data (Figs. 6A and B). The nearest neighbor statistics of all 
three mPIPP mosaics are statistically indistinguishable from each 
other and all three are statistically indistinguishable from the m623 
field (Figs. 6C and D). Figure 6E plots a topological disorder 
parameter [33,40], fi 2 , measuring the spread of the distribution of 
the number of sides in each Voronoi polygon (Fig. 6A, see 
Materials and Methods). \i 2 is plotted for 99 simulated mosaics as 
dots along with the value obtained from the data as a blue 
horizontal line. For each value of y this line falls within the 
distribution of simulated values, indicating a good fit between 
mPIPP model mosaics and the data. Similar results were obtained 
when considering the w81sl mosaic (data not shown). We 
conclude that the local statistics of cell positions is independent 
of the modulation parameter y and matches the statistics of 
experimentally observed mosaics. 

In summary, the proposed modulated PIPP generates semi- 
regular mosaics of ON and OFF center cells with realistic spatial 
statistics that within the statistical connectivity model predict 
experimentally observed OPMs. For increasing y the local spatial 
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Figure 5. Statistical properties of mPIPP mosaics. A mPIPP RGC mosaics obtained with modulation parameter y = 0 (left), y=2.5 (middle) and 
y = 30 (right). ON (OFF) cells are displayed as empty (filled) circles. Dipoles are indicated as colored bars with colors indicating their preferred 
orientation. OPM in the background is the region used to modulate the PIPP. Colorcodes as in Fig. 1C. Parameters for the PIPP were chosen as in [34] 
for cat beta cell mosaic m623 [35]. The cortical magnification is f = 1.7mm,. /mm,.. Dipoles were extracted with tY = 80/<m (see Materials and Methods). 
B Correlation of dipole angles for the mosaics shown in A (red circles: y = 0 (left), 2.5 (middle), 30 (right)). Dark red line indicates correlation function of 
the modulating OPM region in A. C Cross-correlation between dipole orientations and modulating OPM region for different values of the modulation 
parameter-;'. Red line: numerical simulations. Blue line: analytical prediction (see Materials and Methods). D RGC dipole pattern obtained for an mPIPP 
with y = 20 from a 3A x 3A-OPM. All other parameters as in A. For clarity, only the dipoles formed by the mosaics are shown as colored rectangles. 
The modulating OPM region is shown in the background. Black square indicates the OPM region used to simulated the mPIPP in A. E As B for the 
mosaic in D [y = 20). Note, that the dipole correlation function closely follows the correlation function of the modulating OPM. All error bars indicate 
95% confidence intervals of bootstrap distributions. 
doi:1 0.1 371 /journal.pone.00861 39.g005 
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Figure 6. Spatial properties of cell positions in mPIPP mosaics are independent of y and in agreement with experimental data. 

Spatial statistics of the experimental OFF cells in mosaic m623 (column 1) and mosaics obtained by simulating an mPIPP model for three different y 
values (columns 2-4) and fixed £ = 1.7mm t ./mm,. and rf = 80/im. A Central region of the mosaics. Each point denotes one receptive field midpoint and 
is surrounded by its Voronoi polygon. Scale bar: 250 fim. B Autocorrelogram of the points in A, with annuli drawn 25 \.im apart. C Cumulative 
distribution of nearest neighbor distances (G-function, see Materials and Methods). The gray region shows the 95% confidence interval of 
distributions from mPIPP simulations, and the solid line reflects the data (reprinted from column 1). D L-functions for data and models, drawn in the 
same format as for panel C. The dashed line indicates the expectation for Poisson point process (complete spatial randomness), L(f) ~ t. E Topological 
disorder parameter /< 2 for 99 realizations of the mPIPP (black dots) and for the data (horizontal blue line). 
doi:10.1371/journal.pone.0086139.g006 



correlation between dipole orientations increases yet long-range 
positional order remains absent in the generated mosaics (Fig. 5). 
By varying a modulation parameter y, we are able to tune cross- 



correlations between dipole orientations in the RGC mosaics and 
the preferred orientations given in the modulating regions of 
OPMs. The spatial statistics of ON/ OFF cells positions, however, 
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is unaffected by this modulation and all y-values yield mosaics with 
spatial statistics consistent with the data. 

Using mPIPP Mosaics to Assess the Statistical Power of 
Correlation Analysis 

Having characterized its statistical properties, we used the 
mPIPP to estimate the size of a data set necessary to detect 
correlations between dipole angles that could seed the emergence 
of locally smooth and repetitive OPMs. More specifically, we 
wanted to estimate the so-called false negative rate of a test for 
dipole angle correlations, i.e. the rate of rejecting the presence of 
correlations when correlations of a certain degree are present in 
the ensemble the data is drawn from. In principle, testing the 
presence of positive local dipole angle correlations will require 
fewer and smaller mosaics than the detection of presumably 
smaller negative correlations on larger spatial scales. To estimate 
the necessary size of a data set, we generated a number of N 
mPIPP mosaics with known degree of spatial correlation (cf. 
above) and varied their spatial extent. We then statistically 
compared these ensembles to N conventional PIPP mosaics of 
the same spatial extents [34] (see Materials & Methods for details). 
Since PIPP mosaics lack the spatial dipole structure necessary for 
seeding cortical OPMs, this comparison can be employed to 
estimate a lower bound for detectable relevant dipole angle 
correlations, given a data set with N mosaics of fixed size, all 
measured at similar eccentricity. The latter implies that correla- 
tions functions of the mosaics can be averaged to improve the 
statistical power of the test. The two cat beta cell mosaics come 
from different eccentricities. Therefore, their correlation functions 
cannot be averaged. 

Figure 7 depicts the estimates of the false negative error rate fi of 
the statistical test for the presence of dipole angle correlations. 
Figure 7 A shows the false negative rate for the detection of positive 
local correlation in mPIPP mosaics simulated with different 
modulation parameters y as a function of the number of mosaics 
and their area. Even for small y, already one mPIPP mosaic of an 
area of 1 mm 2 (the size of the cat beta cell data sets) yields 
/? = 0.07. This means that 93% of the realizations of such mosaics 
are statistically distinguishable in terms of local dipole angle 
correlations from an ensemble of PIPPs. Similarly, 5 mosaics with 
an area of around 0.4 mm 2 would be sufficient to reliably detect 
even weak local correlations (y = 3). Figure 7B displays the 
probability of failing to detect negative correlation of dipole angles 
around a distance of 300/wz in the mPIPP mosaics. As expected, 
averaging over more mosaics or, alternatively one larger mosaic is 
needed to detect such anti-correlations. From the above analysis, 
we conclude that both cat beta cell mosaics analyzed here are 
sufficiently large to detect even weak local positive correlations of 
dipole angles that could seed iso-orientation domains in VI, if 
present. However, a larger data set (e.g. N>9 mosaics of 1 mm 2 
measured at the same eccentricity) is needed to reliably test for the 
negative correlations that would be indicative of seeding spatially 
repetitive OPMs. Note that, since the mPIPP mosaics generally 
contain slightly fewer dipoles than experimental RGC mosaics (cf. 
Fig. 5A), all of the above estimates of statistical power are 
conservative and most likely even fewer or smaller samples of real 
mosaics would suffice to detect the presence of the respective 
correlations. 

Comparison between mPIPP Mosaics and Cat Beta-cell 
Mosaics 

The absence of spatial correlations in dipole orientations in both 
published cat RGC mosaics can be employed to determine an 



upper bound for the modulation parameter y of the mPIPP. To do 

so, we synthesized mPIPP mosaics for a wide range of cortical 

mm c mm c 

magnification factors c between 0.2 and 2.5 (estimated 

mm, mm,. 

value for both mosaics is {=1.7) and for a wide range of 
modulation parameters y between 0 and 30. The maximum dipole 
distance was fixed at d = 80{im. However, similar results were 
obtained with d = 60 fim and d=\00^m. Each mPIPP mosaic 
realization was modulated by a randomly chosen rectangular 
OPM region from cat area 17. The size of the OPM region was 
determined by the respective magnification factor £ (see Materials 
and Methods). For each pair of values (£, y), we generated 100 
mPIPP realizations and calculated their dipole orientation 
correlation functions. Figures 8A and B depict these correlation 
functions for £ = 1.7mm c /mm r and y = 30 (A) and y = 0 (B). For 
y = 0, dipole orientation correlations of mPIPP mosaics and the 
experimentally measured mosaic perfectly overlap. For y = 30, 
correlation values of experimental data lie far outside the min/ 
max range of correlation values obtained with the mPIPP mosaics. 
To determine the range of y for which mPIPP mosaics are 
consistent with the data, we then asked for each pair of values ({, 
y), how likely the correlation values in the data for distances 
smaller than 200fun are generated by the corresponding mPIPPs. 
For each pair (£, y) this resulted in a Monte-Carlo p-value, 
indicating the likelihood of finding the correlation value of the 
experimental mosaic in the ensemble of mPIPPs. The p-values for 
different (£, y) are depicted in Figures 8C (m623) and D (w81sl). 
For both mosaics, p< 0.05 if y>3, indicating that the experimen- 
tal data is only consistent with very small values of y. For these 
values, positioning of RGC cells in the mPIPP is only very weakly 
modulated by the OPM region (see Figs. 5A and B). In fact, the 
data from both RGC mosaics w81sl and m623 are most consistent 
with y = 0, i.e. the previously described PIPP [34] with complete 
absence of modulation from the OPM region. As emphasized 
before, in this regime the VI OPM predicted by the statistical 
connectivity model does not exhibit a typical column spacing and 
therefore the model cannot account for the experimentally 
observed semiregular structure of visual cortical orientation maps. 

Discussion 

In this paper we tested a fundamental prediction of the 
statistical connectivity model [18-23]. When considered in the so- 
called dipole approximation, i.e. when each VI neurons receives 
input from only a few number of RGCs, the statistical connectivity 
model predicts that RGC dipole angles exhibit spatial correlations 
similar to the spatial correlations of orientation preferences in VI 
OPMs. This means that dipole angles should be locally positively 
correlated and anti-correlated on larger scales. We analyzed two 
cat beta cell mosaics as well as one primate parasol receptive field 
mosaic searching for the presence of such correlations. All three 
mosaics lack spatial correlations on the relevant spatial scales 
(Figs. 2, 3). Weak local correlations can be attributed to receptive 
fields being semi regularly spaced and are not sufficient to seed 
realistic cortical OPMs (Fig. 4). 

To investigate whether the absence of detectable correlations 
was merely a consequence of the small size of the available mosaic 
data sets, we used a novel point process (mPIPP) which generates 
realizations of an ON/ OFF ganglion cell mosaic that could seed 
realistic OPMs [18-20]. By modulating the positioning of ON/ 
OFF cells with experimentally obtained OPMs, the mPIPP 
generates semi-regular ganglion cell mosaics which within the 
statistical connectivity model will generate realistic spatially 
repetitive OPMs. Notably, no Moire-Interference mechanism 
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Figure 7. Estimating the statistical power of the test for the presence of spatial correlations in RGC dipole angles. A False negative 
rate (probability of failing to detect the positive local correlation) for mPIPP mosaics simulated with different modulation parameters y as a function 
of the number of mosaics N and their area size (see Materials & Methods for details). Green box indicates size of cat RGC mosaic data sets analyzed in 
the present study (N = 1, area size ss 1 mm 2 ). Note that even for y = 3, the false negative rate with this data set size is very small. B False negative rate 
for detecting negative correlation of dipole angles mPIPP mosaics around a distance of 300/im (see Fig. 5, see Materials & Methods for details). In all 
panels, a cortical magnification of c = LVmnv/mm,. was assumed and PIPP parameters were taken from the fit to the m623 mosaic [34]. 
doi:10.1371/journal.pone.0086139.g007 



needs to be invoked to generate periodicity of the OPMs. The 
process extends previous work on PIPPs [33,34] to generate ON/ 
OFF cell mosaics with realistic spatial properties. By varying a 
modulation parameter y, the influence of the OPM region on the 
ON/ OFF cells' position could be arbitrarily and predictably tuned 
(Fig. 5). The local spatial statistics of the model mosaics agree well 
with experiment (Fig. 6). They are essentially insensitive to the 
modulation through the OPM. However, the inferred mosaics are 
characterized by a salient local correlation between neighboring 
dipole orientations, anti-correlation at intermediate distances and 
correlation at greater distances (Fig. 5), an effect of the typical 
distance between adjacent columns preferring the same orienta- 
tion in the OPM. 

The mPIPP mosaics were then used as reference mosaics with a 
predictable degree of dipole correlations to determine the 
statistical power of inferring the presence/ absence of correlations 
from finite samples of RGC mosaics. W e find that both cat beta 
cell mosaics are sufficiently large to reliably detect even weak local 
positive correlations of dipole angles that could seed iso-orientation 
domains of the size observed in experimental OPMs (Fig. 7). 
However, a larger data set of mosaics measured at the same 
eccentricity would be needed to reliably detect the anti-correla- 
tions on larger spatial scales that would be indicative of seeding 
spatially repetitive OPMs. These findings then prompted us 



determine an upper bound for the local dipole angle correlations 
in the data (Fig. 8). Our results show that even weak correlations 
are ruled out by the data. Thus, experimentally measured RGC 
mosaics lack the local dipole structure to seed iso-orientation 
domains of the size observed in experimental OPMs. 

As explained above the data sets analyzed in the present study 
are not sufficiently large to rule out a weak periodicity of dipole 
angles in experimentally measured RGC mosaics. However, we 
would like to emphasize that the presence of local correlations in 
ON/ OFF dipole angles appears as a necessary prerequisite for the 
statistical connectivity model to yield spatially repetitive OPMs. 
When hexagonal RGC mosaics are considered, dipole angle 
correlations emerge via Moire Interference of hexagonal mosaics, 
that are positionally ordered over long distance (>lmm) [20,21]. 
Our algorithm to "reverse engineer" plausible realizations of 
RGCs from measured OPMs shows that the statistical connectivity 
model does not necessarily rely on such long-range positional order. 
In the PIPP model, correlations of neighboring ON/ OFF dipole 
angles could instead be build into the mosaics by a simple 
modification of the intra-mosaic interaction function studied in 
[34]. The absence of angular correlations in the data not only 
constrains the modulation parameter y of the mPIPP to very small 
values. It also provides further evidence against the seeding of 
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Figure 8. Constraining the modulation parameter y with experimental data. A Correlation functions of dipole preferred orientations for 100 
mPIPP realizations (pale pink dots, y = 30, cortical magnification £ = 1.7mm ( ./mm,., rf = 80/»«). Black drawn line indicates average correlation function, 
dashed lines show + la deviation from the mean. Red dots indicate correlation function of dipole orientations for the mosaic m623 (redrawn from 
Fig. 2H). Insets show the T-distribution for Monte-Carlo data (blue) to estimate the p-value for the observed value (red) (see Materials and Methods). 
Note the correlation-anti-correlation-correlation structure of the correlation function. B As A but for y=0. C Monte-Carlo p-values (see Materials and 
Methods) for different values of-/ and the cortical magnification factor { for m623. Red line indicates p = 0.05 significance value. D as C but for mosaic 
w81s1. 

doi:1 0.1 371 /journal.pone.00861 39.g008 



cortical OPMs via feed-forward projection of a Moire interference 
pattern from the retina. 

One main simplification of the statistical connectivity approach 
is that the transformation of visual inputs by the LGN is usually 
ignored. RGC afferent terminal axons diverge in the A-lamina, 
providing a one-to-many mapping between one X-RGC axon and 
several X-type geniculate relay cells in cat LGN. Moreover, 
geniculate neurons often receive input from several retinal 
afferents [41]. However, these multiple retinal inputs have mostly 
overlapping receptive-field centers [41]. This has important 
implications for how the pattern of retinal dipoles might be 
transformed into a pattern of LGN dipoles. The increased density 
of LGN dipoles might not interpolate the sparsely sampled dipole 
pattern of the retina. Therefore, the size of the iso-orientation 
domains set by the LGN dipoles does not increase with respect to 
the iso-orientation domains set by the RGCs. Instead domains of 
dipoles of equal orientation sharply confined around RGC dipoles 
(similar to a Voronoi tessellation of the dipole center positions) 
might be the most likely outcome of the thalamic transformation. 
If so, the thalamic transformation would not be able to transform 
the uncorrelated dipole pattern in the retina into a pattern with 
angular correlations on the relevant scales. 

The statistical connectivity model has been advanced as a 
theory explaining the establishment of a blueprint for the VI 
orientation preference map during early visual development based 
on retinal inputs from the contralateral eye [20] . In fact, the spatial 
layout of initially established contralateral eye dominated OPM is 
similar to the binocular OPM of the mature animal [42]. Hence, 



retinal inputs from the contralateral eye alone should be 
sufficiently structured to generate a spatially repetitive OPM with 
typical column spacing in the millimeter range. These findings 
usually justify the analysis of properties of mosaics from one eye 
only while, of course, in the adult most visual cortical neurons are 
binocular [1]. If, as our analysis suggests, mosaics in the individual 
eyes lack the spatial structure necessary to seed realistic OPMs, it 
appears unlikely to us that by combining inputs from ipsilateral 
and contralateral eye, a correctly spatially structured seed could be 
established. To establish the necessary spatial structure from two 
spatially unstructured inputs, a considerable fraction of VI 
neurons would have to receive their OFF subfield input from 
one eye and their ON subfield input from the other. If this were 
the case, the OPMs measured by stimulating the ipsilateral or 
contralateral eye only in the adult animal should be considerably 
different. This is, however, not what Crair et al. have shown [42] . 
We conclude that competing eye inputs are unlikely to offer an 
alternative non-cortical mechanism for generating the structure of 
experimentally measured OPMs. 

The hypothesis (in this study called statistical connectivity 
model) that cortical OPMs could emerge from a spatially 
structured retinal organization has been previously considered 
with three different RGC mosaic classes and in two different 
regimes in terms of the number of feed-forward retinal inputs N 
that a VI neuron samples from. Both, Soodak and Ringach 
considered noisy hexagonal RGC mosaics and a large number of 
inputs to each cortical VI cell (Nx 20) [18,23]. In this setup, the 
statistical wiring model results in OPMs without a typical column 
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spacing and, hence, cannot account for the spatial structure 
observed in experimentally measured OPMs. Park & Ringach 
considered the model with noisy hexagonal RGC mosaics in the 
so-called dipole approximation, where input to a VI cell is 
dominated by on average a close-by pair of one ON and one OFF 
cell [20,21]. In this regime, OPMs predicted by the model exhibit 
a typical spacing between adjacent columns. However, as oudined 
in the introduction, the positional statistics of RGC mosaics is not 
well described by noisy hexagonal lattices [33] and the PIPP 
mosaics provide a much better fit to the data. The statistical 
connectivity model considered with these more realistic mosaics 
fails to generate maps with a typical column spacing [19,33]. Our 
study introduces a third class of RGC mosaics (mPIPP mosaics) 
with realistic positional statistics and a repetitive structure of 
cortical OPMs. The improved match between data and model, 
however, comes at the natural expense of spatial correlations in 
dipole angles in the mPIPP mosaics which are not found in 
experimentally measured mosaics. 

The statistical connectivity model's mismatch with empirical 
evidence in each of the above mentioned parameter regimes 
suggests that the spatial layout of OPMs is not determined by the 
structure of RGC mosaics and may instead result from intracor- 
tical mechanisms. One of the most striking demonstrations of how 
such intracortical mechanisms can shape visual cortex architecture 
comes from cat primary visual cortex. Using pharmacological 
treatments, Hensch and Stryker [43] locally altered the balance 
between intracortical inhibition and excitation during OPM 
formation. Enhancement of inhibitory circuits locally widened 
column spacing in VI while local reduction of inhibition 
broadened the spacing of columns. While those findings are 
difficult to reconcile with the statistical connectivity hypothesis, 
theories in which cortical columns arise from an intracortical 
interplay between inhibition and excitation, e.g. [15,16,44,45], 
could provide simple and plausible explanations for such an effect. 
Similarly, the progressive interareal and interhemispheric match- 
ing of features of columnar architecture such as the local column 
spacing in cat VI suggests a strong influence of activity-dependent 
intracortical and even interareal interactions during postnatal 
column formation and critical period refinement and reorganiza- 
tion [9] (see also [46]). 

The fact that RGC mosaics lack the spatial structure to seed 
realistic cortical OPMs does a priori not rule out the possibility 
that retinal/thalamic receptive field mosaics might exert an 
influence on the layout of OPMs during postnatal development 
or even provide a rough blueprint of visual cortical maps. In fact, 
the representations of retinal blood vessel angioscotomas in the 
visual cortex in some squirrel monkeys are a striking demonstra- 
tion of the influence of retinal features on cortical selectivity 
layouts [47] (see [48] for a theoretical treatment in terms of 
activity-dependent mechanisms). If this is the case, preferred 
orientations of RGC dipoles on the retina might still be a fairly 
good predictor of cortical orientation preference for a substantial 
fraction of neurons, and intracortical network interactions might 
serve to organize a retinal seed of orientation preference into a 
pattern with a typical spacing and a semi regular arrangement of 
pinwheel centers. Furthermore, it remains possible that a higher 
order statistical structure of RGC mosaics beyond simple dipoles 
(or their thalamic transformation) is capable of driving the 
formation of realistic OPMs. Including retinal biases or constraints 
into existing models for the activity-dependent self-organization of 
visual cortical orientation preference is needed to elucidate the 
interplay between subcortical feed-forward constraints and intra- 
cortical network self-organization. 



Finally, our conclusions are necessary limited by the fact that we 
do not have retinal fields and corresponding orientation maps 
from the same animal. This is clearly a large task, and one outside 
the scope of our current paper. Since anatomical RGC mosaics 
were first presented in 1981, there have been only a few mosaics 
published [30,35]. In particular, we have analyzed all retinal 
mosaic data available to us from the two recent publications 
[20,33]. Our results from the theoretical modeling presented here 
should thus also be viewed as predictions that we would like to 
revisit once more suitable data are available. 

Materials and Methods 

Optical Imaging Data & Preprocessing 

Optical imaging data from twelve adult cat area 1 7 hemispheres 
were used in this study. No experiments were carried out for the 
sole purpose of the present study. We reanalyzed partially 
published data, collected by the laboratory of Z. Kisvarday at 
Ruhr University Bochum. The original animal license for these 
experiments was issued to Prof. Ulf Eysel (Dept. Neurophysiology) 
and the research program was supported by SFB509. All 
experiments were conducted according to ethical regulations 
issued by the Ruhr University Bochum and conformed to the 
guidelines of the European Communities Council Directives, 1986 
(86/609/EEC) as well as the German Animal Welfare Act. All 
animals derived from in house animal farm or from registered 
breeders for experimental animals. 

Surgery and preparation protocols have been described in detail 
elsewhere [49,50] . Briefly, optical imaging of intrinsic signals were 
conducted on anesthetized (initial anesthesia: a mixture of 
ketamine, 7 mg/kg Ketanest (Parke-Davis, Berlin, Germany), 
and Xylazine, 1 mg/kg (Rompun, Bayer Belgium, Sint-Truiden, 
Belgium), i.m.; prolonged anesthesia: 0.4-0.6% halothane in a 1:2 
mixture of O2 and N2O using artificial ventilation) and paralyzed 
(alcuronium chloride (0. 15 mg/kg/h, Alloferin, Hoffman-La 
Roche, Grenzach-Whylen, Germany, i.a.) animals using the 
imaging system Imager 2001 (Optical Imaging, Germantown, 
NY) and the data acquisition software VDAO_ (Version No. 
VDAQ218k, Optical Imaging). A craniotomy was performed on 
one hemisphere between stereotaxic coordinates (Horsley-Clarke) 
P7-A12 and L0.5-L6.5 to expose the cortical region correspond- 
ing to the representation of the central and lower parts of the 
visual field in both area 17 and area 18 [27]. Animals were 
monitored continuously throughout all procedures to ensure that 
adequate anesthesia was maintained. Area 17 identification was 
made on the basis of stereotaxic coordinates. Before acquiring 
data, the camera was focused at 650-750/im below the cortical 
surface that was illuminated with 609 + 5 nm light. Visual stimuli 
were presented to one eye. Full-field visual stimuli were presented 
on a video screen (SONY, Pencoed, UK) in 120 Hz noninterlaced 
mode. High contrast, square-wave gratings were generated at 
optimal spatial (0.1-0.2 cycle /deg) and temporal frequencies (1— 
2 Hz), using the stimulus generator systems VSG Series Three 
(Cambridge Research Systems, Rochester, UK). After the 
recording session, the animals were euthanized with an overdose 
of anesthetics (Nembuthal). High vs. low spatial frequency stimuli 
characteristic, respectively, for area 17 and area 18 were used to 
visualize the area border. 

Difference maps were obtained from single-condition maps as 
described previously [49,50]. For each hemisphere a region of 
interest (ROI) was defined containing the imaged part of area 17. 
ROIs typically contained around 10-20 A 2 were A is the typical 
column spacing of the OPM. Each raw difference map was Fermi- 
bandpass filtered as in [10,11]. Low-pass cut-off was chosen as 
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1.5mm, high-pass cut-off was chosen as 0.4mm. This preprocess- 
ing ensured efficient removal of high-frequency noise from the 
CCD-camera and low frequency variations in signal strength while 
only weakly attenuating structures on the relevant scales. 

Experimental Mosaics 

All RGC mosaics used in the present study are available for 
download from the website of Dr. Stephen J. Eglen (http://www. 
damtp.cam.ac.uk/user/ sje30/ data/mosaics/). 

Numerical Procedures for mPIPP Mosaics 

Numerical procedures to obtain mPIPP ON/OFF mosaics were 
chosen as previously described [19,34]. In short, we initially 
positioned «off OFF-cells and «on ON-cells independendy 
according to a two-dimensional Poisson point process with cell 
density matched to the density of cells in m623 and w81sl 
respectively. We then updated these positions according to the 
following loop (see Fig. 9A): For each ON-center cell a new 
candidate position was generated at random. Considering the i-th 
ON-center cell, this new position was accepted with probability 

Pi = 

"ON . "OFF . . (5) 

.^on,on(|x on — x' on |)- n ^on,off(x on ,x' off ), 

where //on.off is defined as in Eq. (3). After updating all ON 
center cells' positions, the procedure was repeated for the OFF 
center cells, here for the cell number i 



of d. Spatial correlations C{R) between dipole angles (j>(x) and 
<^(y) were calculated as [51,52] 

C(K)= (cos(^(x)-0(y))) (J! _ i/2)s||x _ y||2<(J?+A/2) , 

where b is a fixed bin size. To obtain Fig. 5B, Figs. 2G^J and Fig. 8, 
for a given mosaic, the bin size was chosen such that the diagonal 
of the rectangular retinal section considered contained 20 
equidistant bins. Cross-correlation between the modulating 
OPM 9(x) and the dipole orientations <^>(x) (Fig. ID, lower panel) 
were determined with 



C cross = (cos(<Kx)-0(x))) dipoles • 

where the average is taken over all dipoles in the retinal region 
considered. The expected value of C cross for an mPIPP with a 
given y can be estimated using the probability p of positioning an 
ON/ OFF cell 



/<A0)=iexp(y(cos(A0)-l)), 

where AO is the difference between the dipole orientation </) and 
the preferred orientation specified in the orientation map 8(x) and 
A is a normalizing factor, i.e. N= J <i(A0) exp (y( cos (A#) — 1)). 
For large enough retinal regions, we have 



Pi = 

"OFF . "ON . . (6) 

II ./ioff,off(|x off — x off |)- IJ 77 0 n,off(x on ,x off ) 

with different parameters for Aoff.off and Aon.on- Both loops 
were repeated between 20 and 50 times. Note that from one 
complete iteration of the algorithm to the next, the absolute 
positions of all cells show no correlation (cf. Fig. 9). However, the 
spatial statistics of the pattern including dipole angle correlations 
are stable characteristics of the mosaics after only a few iterations. 

Model parameters besides y and d were chosen as to match the 
local spatial statistics of m623, w81sl [34] and g09 [33], 
respectively, and are summarized in Table SI. RGCs close to 
the boundary of the simulated domain can only form dipoles in a 
subset of all possible directions, potentially leading to boundary 
effects in the mPIPP simulations. Since the dipole orientation is 
enforced by a given OPM, we usually observe a slight decrease in 
the dipole density towards the very edge of the retinal patches. 
However, for the rather small beta-values that are consistent with 
the experimental data (cf. Fig. 8), such boundary effects were 
observed to play a negligible role. Source code for generating 
mPIPP mosaics along with Matlab code for visualization is 
provided as supplementary material to the manuscript. 

Dipole Extraction and Correlation Function 

Following [37], we assumed a pair of ON/ OFF ganglion cells at 
position x and y respectively to form a dipole if their distance was 
smaller than a parameter d, i.e. ||x — y|| 2 <d. A dipole's preferred 
orientation was defined as in Eq. (2). Note that this orientation is 
orthogonal to the orientation of the dipole vector connecting the 
ON/OFF pair (see Fig. 1A). A dipole's position was defined as 
(x + y)/2, i.e. halfway between the ON and the OFF cell. Note 
that one RGC can form multiple dipoles depending on the choice 



(cos (A6) ) dipoles = |" ^ d(A8)p(A6) cos (A6) = g| , 

where Ij(y) is the j-th modified Bessel function of the first kind. 
This function is depicted in Fig. 5C (blue line) along with the 
numerically obtained cross-correlation values (red line). 

Cortical Magnification Factor at the Position of the 
Measured Mosaics 

To determine the VI OPM region sizes to which the two 
published cat beta cell mosaics correspond, we estimated the 
cortical magnification factor at the respective retinal positions. The 
center of mosaic w8 1 s 1 was located 1 9 below the mid line of the 
visual streak, 4mm from the area centralis [35]. With respect to a 
cartesian coordinate system with origin at the area centralis and x- 
axis along the visual streak this position is 

4mm,. cos (19 )\ _ / 1.3mm r \ 
4mm r sin(19 ) J \ 3 - 8mm r/ 

In the cat retina, 1mm, roughly corresponds to 4.4 deg visual 
angle [53,54]. Thus 

/ 1.3mm, \ / 5.9° \ 
\3.8mm r J^yn. 2° )' 

From [27], we determined the cortical magnification for 
elevation 5.9 and azimuth 17.2° to be about 0.15mm^/deg 2 . 
Thus, via 
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for i = l...noN 

delete position Xq N 

repeat until new x'q N accepted 

propose new position Xq S 
compute pi according to Eq. (5) 
accept new position w/ probability pj. 



for j = 1...UOFF 

delete position Xq ff 

repeat until new x^ FF accepted 

propose new position Xq ff 
compute pj according to Eq. (6) 
accept new position w/ probability pj 
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Figure 9. Constructing RGC mosaics from measured OPMs with an mPIPP. Upper panel: OPM measured in cat VI. Below, left: Inset of OPM 
used for modulating the mPIPP; right: Initial condition of the Monte-Carlo optimization procedure. ON/OFF RGC mosaic are specified by a 
homogeneous Poisson process. Positioning of RGCs is irregular. The preferred orientations of dipoles (colored bars) do not match the orientation 
preferences of the OPM region. Lower most panel: ON/OFF RGC mosaic after 20 iterations of the Monte-Carlo procedure with y = 30 and 
£ = 0.6mm t ./mm,. for mosaic m623 (see Materials and Methods). Positioning of RGCs is semiregular and preferred orientation of dipoles are almost 
perfectly aligned with orientation preferences of the OPM region (background). Colorcode as in Fig. 1C. B Orientation map region, superimposed 
mPIPP mosaics (left) and corresponding dipole correlation function (right) for 0,5,40, and 80 iterations (upper most to lower most panels). Note that 
the dipole angle correlations settle around their final values within only a few iterations of the process, whereas the precise positions of the ON/OFF 
cells continues to change. 
doi:1 0.1 371 /journal.pone.00861 39.g009 



0.15 



mm? 4 4 deg 
deg 2 mm,. 



■ 1.7- 



mm,- 



we estimate that, at the location of w81sl, 0.6mm on the retina 
correspond to 1mm visual cortex. The typical column spacing of 



cat OPMs is A c «lmm c [39,55]. Hence, we would expect a 
periodicity of about A,. = 0.6mm,. on the retina. 

The mosaic m623 has been obtained from 5 mm eccentricity 
and 5.5 deg below the mid line of the visual streak. This 
corresponds to the point (2.1 ,22.6 ) in the visual field, which 
again has 0.1 5mm 2 /deg 2 cortical magnification. Thus 
0.6mm, = 1mm,. for m623. 
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The mosaic g09 was measured at 9 mm eccentricity (temporal 
retina) at 41 visual angle in Macaca mulatta, Macaca fascicularis. 
At this point 1.4mm r = lmm c ([20], Suppl. Inf.). 

Spatial Statistics of mPIPP Mosaics 

We have analyzed the spatial statistics of mPIPP mosaics using 
previous methods [33,34]. In short, to obtain the G-function, for 
each point in a mosaic the distance to its nearest neighbor was 
calculated and the cumulative distribution of these distances was 
computed. The L-function is the scaled expectation of the number 
of points observed within a given distance of any point. In 
Figures 6C and D, we have drawn the 95% confidence levels from 
99 mPIPP simulations (gray shading). Informally, if the measure 
from the observed mosaics (solid lines in Fig. 6) falls within the 
confidence intervals, then the model is a good fit to the observed 
mosaic. Topological disorder, fi 2 was quantified using [40] 

^2=I]("-6) 2 P„, 

n 

where P„ is the probability of a Voronoi polygon having n edges. 

Calculation of p-values for Local Spatial Statistics 

All p-values used in the present study are Monte Carlo P values 
[34,56]. Each model is run, using the same parameters but with 
different initial conditions, 99 times. Each mosaic (both observed 
and simulated) is then compared against the other 99 mosaics, and 
a certain measure T is calculated. The 99 T values are sorted, and 
the rank of the T measure corresponding to the observed field is 
divided by 99 to calculate a p-value. p-values of 0.05 and smaller 
indicate that the model does not fit the observed data at the 5% 
significance level. The smallest p-value calculated is thus 0.01 with 
99 simulations. 

Calculation of p-values for Spatial Correlation of Dipole 
Angles 

From ./V = 1 00 independent realizations of the mPIPP, we 
obtained circular correlation functions C'(R) with fe[l,100] and R 
being the binned distance (see above). We introduced a measure 

200/otj 

T,= J2 {C'(R)- (C(R))) 2 (7) 

R = 0fim 

where (C(R)) = — y~\ C'(R) is is the average correlation of 

the mPIPP simulations. Note that only bins with R<200fim were 
used (the region of strong positive correlation). The distribution of 
values T quantifies how much the test ensemble of C; deviates 
from the average. Next, we repeated this evaluation with the 
correlation function M(R) obtained from the cat mosaics: 



S= £ (M(R)-(C(R))) 2 (8) 

_R = 0/mj 

From the distribution of the Tj's and value of S, the Monte- 
Carlo p-value was estimated as described above. 

Estimation of the Statistical Power of Hypothesis Test 

The null hypothesis of the statistical test is that RGC dipole 
angles are correlated in space. To estimate the statistical power 
of our hypothesis test for the presence of correlations, it is 



necessary to estimate the probability that the test will reject the 
presence of correlations when the alternative hypothesis is true, 
i.e. the RGC mosaics stems from an ensemble of correlations. 
To do so, we simulated an ensemble of mPIPP mosaics with 
known correlational statistics. Different realizations of the 
mPIPP were then compared to a PIPP control ensemble with 
respect to their dipole angle correlations. In some mPIPP 
realizations, the angular correlation of the mPIPP will be within 
the range of correlations of the PIPP ensemble. Comparing 
these correlations values and concluding that the mPIPP 
realization stems from a PIPP ensemble without angular 
correlations would constitute a type II statistical error (false 
negative). The probability of such an error, i.e. false negative 
rate, is a measure of the statistical power of the test. This 
probability was estimated by repeatedly comparing mPIPP 
realizations to PIPP ensembles. 

More precisely, we first chose a value y, an area size A and a 
number of mosaics N. Within the area A of the retinal patch, 
the RGC density was fixed to the values of mosaic m623, 
PoN = 67mm~ 2 and p OFF = 74mm~ 2 . Then, we calculated 500 
realizations of mosaics for this area and y = 0 to obtain the 
PIPP ensemble. From this ensemble, we randomly drew 100 
mosaics with replacement and calculated their correlation 
functions. From the 100 correlation functions, N were drawn 
and averaged, in total a 100 times. From this set of averaged 
correlation functions we calculate the distribution of T( 
according to Eq. (7). Next, we compare these Tj with an S 
(Eq. (8)) obtain from averaging correlations functions of N 
random mPIPP realizations of the chosen y and with the same 
A. This resulted in a single p-value for this particular set of 
realizations. To estimate the probability of falsely rejecting the 
presence of angular dipole correlation, we calculated a 
distribution of such p-values by repeating the above stated 
steps a 1000 times, p-values larger than 0.05 (our significance 
level) in this distribution indicate the occurrence of type II 
statistical error. Fig. 7A shows the percentage of p-values larger 
than 0.05 (i.e. the false negative rate, ft) as function of A and N 
for different values of y. Fig. 7A shows the same analysis for 
periodicity. It is done analogously with correlations evaluated 
between 200/.WZ and 400/xm. 

Supporting Information 

Figure SI Spatial properties of ON cell positions in 
mPIPP mosaics are independent of y and in agreement 
with experimental data. Same as Fig. 6, but for ON cells in 
RGC mosaic m623. 
(TIF) 

Table SI Model parameters for PIPP and mPIPP 
simulations used in the current study. Parameters are 
chosen such that PIPP mosaics match the local spatial statistics of 
m623, w81sl (see [34]) and G09 (see [33]). 
(PDF) 

Source Code SI Source code for simulating mPIPP 
mosaics 

(C) 

Source Code S2 Source code for simulating mPIPP 
mosaics 

(DAT) 

Source Code S3 Source code for simulating mPIPP 
mosaics 

(M) 
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